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ABSTRACT 

We present a carefully selected sub-sample of Swift Long Gamma-ray Bursts (GRBs), that is com- 
plete in redshift. The sample is constructed by considering only bursts with favorable observing 
conditions for ground-based follow-up searches, that are bright in the 15-150 keV Swift/BAT band, 
i.e. with 1-s peak photon fluxes in excess to 2.6 ph s~ 4 cm~ 2 . The sample is composed by 58 bursts, 
52 of them with redshift for a completeness level of 90%, while another two have a redshift constraint, 
reaching a completeness level of 95%. For only three bursts we have no constraint on the redshift. 
The high level of redshift completeness allows us for the first time to constrain the GRB luminosity 
function and its evolution with cosmic times in a unbiased way. We find that strong evolution in 
luminosity (<5; = 2.3 ± 0.6) or in density (8 a = 1.7 ± 0.5) is required in order to account for the obser- 
vations. The derived redshift distribution in the two scenarios are consistent with each other, in spite 
of their different intrinsic redshift distribution. This calls for other indicators to distinguish among 
different evolution models. Complete samples are at the base of any population studies. In future 
works we will use this unique sample of Swift bright GRBs to study the properties of the population 
of long GRBs. 

Subject headings: gamma-ray: burst - stars: formation - cosmology: observations. 



1. INTRODUCTION 

Gamma-ray bursts are powerful flashes of high-energy 
photons occurring at an average rate of a few per day 
throughout the Universe. They are detected at all red- 
shifts, from the local Universe up to the extreme high 
redshifts (Salvaterra et al. 2009; Tanvir et al. 2009; Cuc- 
chiara et al. 2011a). Our knowledge of the distribution 
of long GRBs through cosmic times is still hampered by 
the fact that most of the observed Swift GRBs are with- 
out redshift. Indeed, the measure of the distance has 
been secured for only ~ 1/3 of the cases. Given the low 
completeness level in redshift determination, the effect of 
possible observational biases could be important in shap- 
ing their redshift distribution (Fiore et al. 2007). This 
fact strongly limit the possibility of well grounded sta- 
tistical studies of the rest-frame properties of long GRBs 
and their evolution with cosmic time. Therefore, it is of 
paramount importance to obtain an unbiased complete 
sample of GRBs, capable to fully represent this class of 
object. 

To this end, we present in this paper a well selected 
sub-sample of the full Swift database. We select bursts 
that have favorable observing conditions for redshift de- 
termination from ground and that are bright in the 15- 
150 keV Swift/BAT band. We find 58 bursts matching 
our selection criteria with a completeness level in redshift 
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determination of 90%. The completeness level increases 
up to ~ 95% by considering the redshift constraints im- 
posed by the detection of the afterglow or host galaxy 
in some optical filters. Therefore, our selection criteria 
allow us to construct a sizable sample of long bursts that 
is (almost) complete in redshift, providing the solid basis 
for the study of the long GRB population in an unbiased 
way. In particular, since our selection is based on the 
brightness in the Swift/BAT band, our sample is not bi- 
ased against the detection of dark bursts, thus providing 
a complete description of the whole long GRB popula- 
tion. 

In the present paper, we will take advantage of the high 
completeness level of our sample to constrain the GRB 
luminosity function (LF) and its evolution with cosmic 
time. In the last few years, this problem has been faced 
by many different authors ((e.g. Porciani & Madau 2001, 
Firmani et al. 2004, Guetta et al. 2005, Natarajan et al. 
2005, Daigne et al. 2006, Salvaterra & Chincarini 2007, 
Salvaterra et al. 2009b, Butler et al. 2010, Wanderman 
& Piran 2010, Campisi et al. 2010, Qin et al. 2010, 
Virgili et al. 2011, Robertson & Ellis 2012). There is a 
general agreement about the fact that GRBs must have 
experienced some sort of evolution through cosmic time, 
whereas the nature and the level of such evolution is still 
matter of debate. Most of the previous works relied on 
the assumption that bursts lacking of redshift measure- 
ments follow closely the redshift distribution of bursts 
with known z. In the past, we tried to overcome this 
assumption by deriving conservative lower limit on the 
level of evolution on the basis of the number of bursts 
detected at z > 2.5 (Salvaterra & Chincarini 2007) and 
of bursts with peak luminosity L > 10 53 erg s _1 (Sal- 
vaterra et al. 2009b). For the first time, thanks to our 
well selected, complete sub-sample of Swift GRBs, we 
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can tackle this issue in a unbiased way. In future works 
we will use this sample to study the correlation between 
physical parameters of the bursts, the properties of the 
burst light curves and of the environment in which they 
explode. 

This paper is organized as follow. In Sect. 2 we de- 
scribe our selection criteria and present our sample. We 
present our models of the GRB LF and redshift distribu- 
tion in Sect. 3, while the results are given in Sect. 4. We 
extrapolate our findings to the detection limit of Swift in 
Sect. 5. Finally, in Sect. 6 we draw our conclusions. 

2. THE SAMPLE 

About 1/3 of all GRBs observed by the Swift satellite 
(Gehrels et al. 2004) has a measured redshift. While this 
represents an enormous improvement with respect to the 
pre- Swift situation, the sample is still far to be consid- 
ered as complete. Jakobsson et al. (2006) proposed a 
series of criteria in order to carefully select long GRBs 
which have observing conditions favorable for redshift 
determination. In particular, they required that: i) the 
burst has been well localized by Swift/XRT and its co- 
ordinate quickly distributed; ii) the Galactic extinction 
in the burst direction is low {Ay < 0.5); hi) the GRB 
declination is -70° < <5 < 70°; iv) the Sun-to-field dis- 
tance is #sun > 55°; v) no nearby bright stars is present. 
While none of the above criteria is expected to alter sig- 
nificantly the redshift distribution of observed GRBs, the 
completness level is increased to ~ 53%0. Still this level 
of completeness is way too low to permit robust popula- 
tion studies. 

In order to construct a more complete sample we re- 
strict ourself to GRBs that are relatively bright in the 
15-150 keV Swift/BAT band. In particular, we select 
bursts matching the above criteria and having 1-s peak 
photon flux P > 2.6 ph s _1 cm~ 2 . This corresponds to 
an instrument that is ~ 6 times less sensitive than Swift. 
58 GRBs match our selection criteria and are listed in 
Table 1 up to May 2011. 52 of them have measured red- 
shift so that our completeness level is 90%. Of these 52, 
all but two (namely GRB 070521; Perley et al. 2009 and 
GRB 080602; Rossi et al. 2012) have spectroscopic con- 
firmed redshift cither from absorption lines over-imposed 
on the GRB optical afterglow or from emission lines of 
the GRB host galaxy. Moreover, for 3 of the 6 bursts 
lacking measured z the afterglow or the host galaxy have 
been detected in at least one optical filter, so that ~ 95% 
of the bursts in our sample have a constrained redshift. 
We note that, while our sample represents only ~ 10% of 
the full Swift sample, it contains more than 30% of long 
GRBs with known redshift. 

The redshift distribution of the bursts in our sample is 
shown in Fig. [TJ In spite of the rather severe cut in the 
observed photon flux, the bursts in our sample have a 
broad distribution in redshift. The mean (median) red- 
shift of the sample is 1.84±0.16 (1.64±0.10) with a long 
tail at high-z extending, at least, up to z = 5.47. 

3. MODEL DESCRIPTION 

The expected redshift distribution of GRBs can be 
computed once the GRB LF and the GRB formation 

6 Up to May 2011, the sample consists in 248 
long GRBs. 132 with measured red shift. See 
http: / / www.raunvis.hi.is/~pja/GRBsample. html 



history has been specified. We briefly recap here the 
adopted formalism and refer the interested reader to Sal- 
vaterra & Chincarini (2007) and Salvaterra et al. (2009b) 
for the model details. 

The observed peak photon flux, P, in the energy band 
Emm < E < E max , emitted by an isotropically radiating 
source at redshift z is 



(l + z)J^l—S(E)dE 
Andl(z) 



P = 



(1) 



where S(E) is the differential rest-frame photon luminos- 
ity of the source, and di,(z) is the luminosity distance. 
To describe the typical burst spectrum we adopt a Band 
function with low- and high-energy spectral index -1 and 
-2.25, respectively (Band et al. 1993; Preece et al. 2000; 
Kaneko et al. 2006). The spectrum normalisation is ob- 
tained by imposing that the isotropic-equivalent peak lu- 
minosity is L = j/ugy eV ES{E)dE. In order to broadly 
estimate the peak energy of the spectrum, E p , for a given 
L, we assumed the validity of the correlation between E p 
and L (Yonetoku et al. 2004; Ghirlanda et al. 2005, 
Nava et al. 2011). 

Given a normalised GRB LF, 4>(L), the observed rate 
of bursts with peak flux between P\ and Pi is 



dt Jo dz 47r 1 + z 



»L(P 2 ,z) 

< / dL'<f>(L'), 

IL{P u z) 



(2) 



where dV(z)/dz — 4iTcd 2 L (z)/[H(z)(l + z) 2 ] is the comov- 
ing volume elemenlQ, and H(z) = Ho [Om(1 + z) 3 + Q\ + 
(1 — f^M — ^a)(1 + - z ) 2 ] 1/ ' 2 - Afi s is the solid angle covered 
on the sky by the survey, and the factor (1 + z) _1 ac- 
counts for cosmological time dilation. Finally, 'I'grb(z) 
is the comoving burst formation rate. 

We explore two general expression for the GRB LF: 
a single power-law with an exponential cut-off at low 
luminosity (exponential LF) and a broken power-law LF. 
The former is described by: 



<j)(L) cx 



and the latter by: 



L 



exp 
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<f>(L) 



( Lcut ) 
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for L < Pr- 



ior L > L c 



(4) 



where L cut is the cut-off (break) luminosity and £j, and 
£/ are the bright- and faint-end power-law index, respec- 
tively. The GRB LF are normalized to unity. In order 
to prevent the integral to diverge we adopt a minimum 
GRB luminosity L m ; n = 10 49 erg s" 1 . We test that our 
results do not change if minimum luminosity of 10 48 erg 



7 We adopted the 'concordance' model values for the cosmolog- 
ical parameters: h = 0.7, f! m = 0.3, and SIa = 0.7. 
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Table 1 

List of the bursts matching our selection criteria. Redshifts or limits are provided in the following references: [1] Berger et al. 2005, [2] 
Fynbo et al. 2009a and references therein, [3] Perley et al. 2009, [4] Foley et al. 2005, [5] Levan et al. 2007, [6] Cucchiara & Fox 2008, [7] 
Rossi et al. 2012, [8] Berger et al. 2008, [9] Berger & Rauch 2008, [10] Kuin et al. 2009, [11] Cucchiara et al. 2008, [12] de Ugarte Postigo 
et al. 2009a, [13] D'Avanzo et al. 2009, [14] Chornock et al. 2009, [15] Perley et al. in prep., [16] Wiersema et al. 2009a, [17] de Ugarte 
Postigo et al. 2009b, [18] Fynbo et al. 2009b, [19] Chen et al. 2009, [20] Xu et al. 2009, [21] Cucchiara et al. 2009, [22] Wiersema et al. 
2009b, [23] Milvang-Jensen et al. 2010, [24] Flores et al. 2010, [25] Cucchiara et al. 2011b, [26] de Ugarte Postigo et al. 2011, [27] 

D'Avanzo et al. 2011. 



(a) based on VLT/X-shooter spectra of the host galaxies obtained within the program 087.A-0451 (PI: H. Flores). The spectra 
were reduced using the X-shooter data reduction pipeline version 1.3.7 (see Goldoni et al. 2006). In the NIR arm spectrum of 
GRB 060306 we identified at the afterglow position the [Oil] doublet emission at z = 3.5, while in the NIR arm spectrum of 
GRB 081221 the [OIII] and Ha emission lines are present at z = 2.26 at the afterglow position. 

(b) for this GRB a redshift of z = 0.84 was reported by Thoene et al, (2007). Images and spectra at the afterglow position 
have been taken also using the VLT/FORS within the program 177.A-0591 (PI: J. Hjorth). We downloaded these data from 
the ESO Archive and reduced them with standard procedures using the ESO-MIDAS package. At the afterglow position we 
could identify an object, showing a continuum signature in the spectra, that we therefore consider to be the host galaxy of 
GRB 060814. We can also identify the galaxy reported by Thoene et al. (2007), at z — 0.84, but this object is offset from the 
afterglow position. VLT/X-shooter spectroscopy of the host galaxy has been performed within the program 084.A-0303 (PI: J. 
Fynbo). We reduced these spectra using the X-shooter data reduction pipeline version 1.2.0 (see Goldoni et al. 2006). Thanks 
to the identification in the NIR arm of the [Oil] doublet, [OIII] A5007 and Ha emission lines associated with the host galaxy, 
we can establish a redshift of z = 1.92 for GRB 060814. 

(c) based on the R— band host galaxy detection in ESO-VLT/FORS2 imaging data obtained within the program 177.A-0591 
(PI: J. Hjorth), taken from the ESO Archive. 

(d) photometric redshift on the bases of the most probable host galaxy association in the XRT error circle (Rossi et al. 2012). 



s 1 is adopted (apart from the value of the normalization 
Vo)- 

4. MODEL RESULTS 

We optimize the value of the model free parameters, 
that is the GRB LF, the normalization rjo and the evo- 
lution parameter, by minimizing the C-stat function 
(Cash 1979) jointly fitting the observed differential num- 
ber counts in the 50-300 keV band of BATSE (Stern et 
al. 2001, 2002) and the observed redshift distribution of 
bursts in our sample with photon fluxes in excess to 2.6 
ph s- 1 cm" 2 in the Swift 15-150 keV banrfl For BATSE, 
we adopt 9.1 yrs of observation with an average exposure 
factor of 0.47, including both Earth-blocking and aver- 
age duty cycle for useful 1.024s continuos record (Stern 
et al. 2002). While our complete Swift sample provides a 
powerful test for the existence and the level of evolution 
of the long GRB population with redshift, the fit to the 
BATSE number counts allows us to obtain the present 

8 For those GRBs lacking of redshift measurement, we randomly 
assign a redshift from a flat z-distribution (taking into account the 
available redshift constraints) not to introduce any a-priori bias. 
We produced 1000 of such random realizations. 



day GRB rate density and to better constrain the GRB 
LF free parameters. It is worth to note that the best-fit 
parameters provide a good fit also of the Swift differen- 
tial peak-flux number counts once the 15-150 keV band, 
the FOV of 1.4sr and the observing lifetime of Swift are 
considered (see also Salvaterra & Chincarini 2007). The 
best fit parameter values together with their l-cr confi- 
dence leve0 for different models are provided in Table 2. 
The corresponding redshift distributions for bursts with 
P > 2.6 ph s^ 1 cm^ 2 are shown in Fig. 1. We test for 
each model the two different GRB LF parametrizations 
described in the previous section and we report in Fig. 1 
the one that gives the best result. 

4.1. No- evolution model 

In a first simple (no-evolution) scenario we can assume 
that long GRBs traces the cosmic star formation and 
that their LF is constant in redshift (L cut (z) = L cnt ,o)- 
In this case the cosmic GRB formation rate is 'I'grb^) = 

9 The errors at l-cr confidence level on the parameters of interest 
(GRB LF and evolution parameter) adopting a C-stat increment 
of 2.30, 3.53 and 4.72 for 2, 3, and 4 parameters of interest, respec- 
tively. 
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Table 2 

Best fit parameters for different models. Errors show the 1-cr confidence level for the parameters of interest (see text in section 4 for the 
details). The last two column report the total C-stat value (i.e. the sum of the C-stat values obtained from the fit of the BATSE and 
Swift dataset) and the Akaike Information Criterion (AIC) score, respectively. We note that in order to properly compare different 
models the AIC criterion has to be considered, where oxp((AJC„,j n — AICj)/2) can be interpreted as the relative probability that the i-th 
model minimizes the (estimated) information loss with respect to the model with the minimum AIC, AIC m i n . The total number of data 



points in the fit is 33. The GRB formation rate at z 



0, r]o, is given in units of Gpc 
units of 10 51 erg s . 



yr and the characteristic luminosity to, 51 is 



rjo^i,(z), where is the normalized cosmic star for- 

mation rate (SFR) and t/q is the present-day GRB for- 
mation rate density in units of Gpc~ 3 yr . We adopt 
the cosmic SFR recently computed by Li (2008) that ex- 
tended the previous work by Hopkins & Beacom (2006) 
to higher redshifts. 

This no-evolution scenario (dashed line in Fig. [1} 
clearly does not provide a good representation of the 
observed redshift distribution of our sample, confirming 
previous findings (e.g. Daigne et al. 2006, Salvaterra & 
Chincarini 2007, Salvaterra et al. 2009b, Qin et al. 2010, 
Wanderman & Piran 2010, Virgili et al. 2011). In par- 
ticular, the peak of the GRB redshift distribution is at 
lower redshift than observed and consequently the rate 
of GRBs at high-z is underpredicted. This is confirmed 
by a more detailed statistical analysis. Indeed, on the 
basis of the Akaike information criterion (Akaike 1974) 
we can safely discard this model being ~ 10 -14 times as 
probable as the luminosity evolution model to minimize 
the information loss (the density evolution model with 
the broken power-law LF is 0.24 as probable as the lu- 
minosity one with the cut-off LF). Moreover, a KS test 
between the no evolution best-fit model and the data of 
our sample gives a chanche probability of ~ 5 x 10~ 5 
that the two distribution are drawn from the same par- 
ent population. 

In the following sections, we will consider evolution 
scenarios that may enhance the number of detections at 
high-z. In particular, we explore: i) a luminosity evolu- 
tion model in which high-z GRB are typically brighter 
than low-z bursts and, ii) two density evolution models, 
both leading to an enhancement of the GRB formation 
rate with redshift. Hybrid models, with both luminosity 
and density evolution, are in principle possible. However, 
the fit with hybrid models results to be very degener- 
ate and does not provide usefull constraints. Therefore, 
we prefer here to consider the two scenarios separately 
to highlight possible similarities/differences between the 
two kinds of evolution. 

4.2. Luminosity evolution model 

Evolution in the GRB LF can provide an enhancement 
of the high-z GRB detection, representing a viable way 
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Figure 1. Normalized redshift distribution of GRBs with P > 
2.6 ph s _1 cm" 2 . Data points (red) show the observed redshift 
distribution and the error bars show the Poisson uncertainties on 
the number of detection in the redshift bin. The dashed line (blue) 
shows the expected distribution for the no-evolution case. Results 
of luminosity and density evolution models are shown with the 
light blue and dark shaded orange regions, respectively, taking into 
account the errors on the evolution parameter. (A color version of 
this figure is available in the online journal.) 

to reconcile model results with the observations. Here, 
we consider the possibility that the cut-off (break) lumi- 
nosity is an increasing function of the redshift, that is 
L cu t{z) — £ C ut,o(l + z) S ' ■ We find that a strong luminos- 
ity evolution with Si — 2.3 ± 0.6 is required to reproduce 
the observed redshift distribution of the bursts in our 
complete sample (light shaded area in Fig. [1} . The re- 
sult does not depend on the assumed expression of the 
GRB LF. 

4.3. Density evolution models 
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Figure 2. Intrinsic redshift distribution of long GRBs for differ- 
ent density evolution models. Dark blue and light yellow shaded 
areas show the results for the density evolution model and for the 
metallicity threshold model, respectively. Dotted (dashed) lines 
reports the modified SFR with Z\ = 2.5 (02 = 2.4). In all cases, 
no evolution of the GRB LF has been assumed. (A color version 
of this figure is available in the online journal.) 

An increase of the rate of GRB formation with redshift 
(on the top of the known evolution of the SFR density) 
will also lead to an enhanced detection of bursts at high- 
z. As a general case, we parametrize the evolution in 
the GRB formation rate as rj{z) = 770(1 + z) Sn . By fit- 
ting our datasets we find that strong density evolution is 
required with 5 n = 1.7 ± 0.5. The amount of evolution 
does not depend on the assumed expression of the GBR 
LF. However, we note that the cut-off LF tends to un- 
derestimate the number of low-z bursts with respect to 
the observed one, leading to some discrepancy with the 
first data point. 

The large value of S n implies an important shift of 
the peak of the GRB formation rate towards higher red- 
shifts with respect to stars. We further investigate this 
issue by applying a "correction" to the shape of the cos- 
mic SFR. This is usually parametrized as three power- 
laws (Hopkins & Beacom 2006, Li 2008) with power- 
index oti = 3.3, 012 — 0.055, ui = —4.46 and breaks 
at zi = 0.993 and z 2 = 3.8 (Li 2008). We fit our datasets 
by letting one of the above parameters free to vary in ad- 
dition to those describing the GRB LF. We find that the 
observed redshift distribution of bursts in our sample can 
be explained either by an increase of the redshift of the 
first break to z\ = 2.5 ± 0.5 or by a hardening of the sec- 
ond power-law to ct2 = 2.4 ±0.4. In both cases, the GRB 
formation rate is found to peak at a much higher redshift 
with respect to stars. The intrinsic redshift distribution 
of GRBs is shown in Fig. [5J 

4.4. Metallicity threshold models 

A subclass of density evolution models foresee the for- 
mation of long GRBs preferentially in low-metallicity 
environments. In this case GRBs will be biased trac- 



ers of the star formation activity being their formation 
suppressed at low- z where most of the galaxies are rel- 
atively metal- rich. Following Langer & Norman (2006), 
we model the fractional mass density belonging to metal- 
licity below a given threshold, Z t h as 



£(*) 



r(O.84,(Z th /Z ) 2 lO°^) 
T(0.84) 



(5) 



where F (r) are the incomplete (complete) gamma func- 
tion, and T(0.84) ~ 1.122. The GRB formation rate is 
then given by *grb(^) oc S(z)**(z). 

We fit our datasets letting Z t h free to vary. The 
available data are well described by models with Z t h < 
0.3 Zq, almost independently on the assumed LF. The 
resulting LF is similar to the one obtained for density 
evolution model. Indeed, the two predicted redshift dis- 
tributions match each others within the uncertainties and 
we refer to the dark shaded curves in Fig. [1] also for the 
metallicity threshold model. 

The range of values for Z t h found in our analysis is in 
agreement with the expectation of the collapsar model 
(Woosley & Heger 2006; Fryer et al. 1999). However, 
such strong metallicity cut-offs seem to be inconsistent 
with the observed properties of GRB host at z < 1 (Man- 
nucci et al. 2011; Kocevski et al. 2011; Campisi et 
at. 2011a). In particular, Campisi et al. (2011a) have 
shown that in the presence of a strong metallicity cut- 
off for the GRB progenitor star, the expected distribu- 
tion of GRB host galaxies in the M-Z and Fundamental 
Metallicity Relation planes is much flatter than observed. 
Larger metallicity thresholds will require some luminos- 
ity evolution in order to reproduce the available data. 
For Z t h = 0.5 Zq, the typical burst luminosity should 
increase with redshift as (1 + z) 1 ' 3±0 - 6 . 

5. SWIFT REDSHIFT DISTRIBUTION 

We compute the redshift distribution expected for the 
full Swift dataset assuming a photon flux limit of P = 0.4 
ph s _1 cm~ 2 and fixing the model free parameters to the 
values given in Table 2. The results are shown in Fig. [3] 
for the different evolution scenarios here explored. The 
models are compared with the redshift distribution in- 
ferred from the sample of GRBs observed by GROND 
(Greiner et al. 2011). This sample has a completeness 
level similar to our but is smaller in size and cover a 
broader redshift range. We find that our evolution mod- 
els provide a good description of the observed redshift 
distribution of the GROND sample without the need of 
any adjustment of the free paramters (a KS test gives a 
probability of 50%), whereas the no evolution model is 
excluded (probability of 5 x 10 -4 ). This further confirms 
the reliability of our analysis strenghtening our conclu- 
sions. 

The predicted redshift distribution of bursts detectable 
by Swift presents a steep raise at low-z peaking at z ~ 2 
with a tail extending at higher redshifts. The median 
redshift of the distribution is z — 2.05 ± 0.15 where the 
error takes into account the uncertainties on the evolu- 
tion parameters. We predict that 3 — 5% of the bursts lie 
at z > 5, consistently with the observational estimate of 
5.5 ± 2.8% (Greiner et al. 2011). This further confirms 
that the majority of dark GRBs are not high-z sources 
but more likely obscured by dust (Perley et al. 2009, 
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Figure 3. Normalized rodshift distribution of GRBs detectable 
with Swift, i.e. with observed photon flux in excess to Pn m = 0.4 
ph s~ x cm -2 . The data show the observed redshift distribution 
as reported by Greiner et al. (2011), on the basis of the almost 
complete sample of GRBs detected by GROND. Error bars show 
the Poisson uncertainties on the number of objects in the redshift 
bin. Model results are shown as in Fig. [T] No attempt to fit the 
observed redshift distribution has been done. (A color version of 
this figure is available in the online journal.) 

Greiner et al. 2011). At z > 8, we expect 1.3 — 3.5 
GRBf3 among the 530 Swift GRBs. This is consistent 
with the two detections reported so far, i.e. GRB 090423 
at z = 8.2 (Salvaterra et al. 2009; Tanvir et al. 2009) 
and GRB 090429B at z ~ 9.4 (Cucchiara et al. 2011). 

It is worth to note that the two evolution scenarios here 
explored predict very similar distributions. Therefore we 
can not distinguish between luminosity and density evo- 
lution simply on the basis of the Swift observed redshift 
distribution. 

6. CONCLUSIONS 

We select a sub-sample of Swift long GRBs that is com- 
plete in redshift. The sample is composed by bursts with 
favourable observing conditions and with 1-s peak pho- 
ton fluxes P > 2.6 ph s _1 cm~ 2 . It contains 58 bursts 
with a completeness level of ~ 90% and provides the basis 
for statistical studies of the properties of long GRBs and 
their evolution with redshift in a unbiased way. GRBs 
can be used to study fundamental issues in astronomy 
and astrophysics, such as the star formation rate and 
the stellar and metal abundances evolution. They can be 
use as tracers of the galaxy evolution, of the ISM com- 
position and to investigate the early universe. Complete 
and fully representative samples of GRBs are therefore 
unique tools to perform these investigations. 

Here, we use the observed burst redshift distribution 
of our complete sample to probe and constrain the evolu- 
tion of the long GRB population in redshift. We confirm 

10 We do not consider here the possible contribution of PopIII 
GRBs that may provide additional GRBs at very high-z (Campisi 
et al. 2011b; de Souza et al. 2011). 



that GRBs must have experienced some sort of evolu- 
tion being more luminous or more numerous in the past 
than observed today. We found that in order to match 
the observed distribution, the typical burst luminosity 
should increase as (1 + z ) 2 - 3±0 - 6 or the GRB rate density 
as (1 + 2;) 1 - 7±0 - 5 on the top of the known cosmic evo- 
lution of the SFR. This result does not depend on the 
assumed expression of the GRB luminosity function. We 
also explore models in which GRBs form preferentially in 
low-metallicity environments. We find that the metallic- 
ity threshold for GRB formation should be lower than 0.3 
Zq in order to account for the observations assuming no 
evolution of the GRB LF. This value, while consistent 
with the expectations of collapsar models, seems to be 
at odd with the observed properties of z < 1 GRB hosts 
(Levesque et al. 2010; Mannucci et al. 2011; Campisi et 
al 2011a). 

Extrapolating our results to P = 0.4 ph s _1 cm~ 2 , 
the predictions of evolution models are consistent with 
available observational constrain without the need of any 
adjustment of the model free parameters. We also note 
that the predicted redshift distributions at the sensitivity 
of Swift for different evolution models are very similar, 
implying that it would be extremely difficult to distin- 
guish among different kind of evolutions only on the basis 
of the observed redshift distribution. 

On the bases of our models we predict that 3—5% of the 
bursts detected by Swift lie at z > 5 consistently with the 
most recent observational estimate of 5.5 ±2.8% (Greiner 
et al. 2011). This indicates that high-z bursts can con- 
tribute only marginally to the observed fraction of dark 
bursts. Finally, we expect 30-50 bursts per year over the 
entire sky with luminosities in excess to 10 53 erg s _1 and 
2-5 with luminosities exceeding 10 54 erg s" 1 . Extreme 
luminous GRBs, such as the current Swift record holder 
GRB 080607 with L iso = 2.26 x 10 54 erg s" 1 (Perley et al. 
2011), should be rare, exploding once every 0.5-1 years 
in the Universe. This corresponds to one detection every 
~ 0.7 — 1.5 years in the field of view of Fermi/GBM (and 
one every 5-12 years for Swift). 
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